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ABSTRACT 

We study the effect of the neutron star spin ~ kick velocity ahgnment observed in 
young radio pulsars on the coalescence rate of binary neutron stars. Two scenarios of 
the neutron star formation are considered: when the kick is always present and when 
it is small or absent if a neutron star is formed in a binary system due to electron- 
capture degenerate core collapse. The effect is shown to be especially strong for large 
kick amplitudes and tight alignments, reducing the expected galactic rate of binary 
neutron star coalescences compared to calculations with randomly directed kicks. The 
spin- kick correlation also leads to a much narrower NS spin-orbit misalignment. 
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1 INTRODUCTION 

There is an increasing interest of a broad astrophysical 
community to coalescing binary compact stars as primary 
sources of gravitational waves for the ground-based grav- 
itational wave observatories. Double neutron stars (DNS), 
observed as binary pulsars, remain to be the most reliable 
objects for gravitational wave searches. On-going LIGO sci- 
ence runs l| Abbot et al.ll2007l ) have already set first exper- 
imental upper limits on their galactic rates of a few per 
year. Astrophysical estimates of the DNS coalescence rate, 
which are based on the binary pulsar statistics or can be 
obtained from population synthesis simulations, are model- 
dependent and vary within more than an order of magni- 
tude around the value 10~^ per year (see rece nt reviews 
iPostnov fc Yungelson|[200^ : iKalogera et aill2007l and refer- 
ences therein). 

The kick velocity imparted to a newborn neutron star 
is an important phenomenological parameter of the core col- 
lapse supernovae and represents one of the major uncertain- 
ties in the theory of binary star evolution. The origin of 
the kicks remains unclear and a number of physical models 
have been suggested (see, for example, [Lai 2004 and ref- 
erences therein). For post-supernova evolution of a binary, 
both the amplitude of the kick and its space direction are 
important. The distribution of the kick amplitudes is usu- 
ally o btained from the analysis of radio pulsar proper mo- 
tions (|Hobbs et al., 2005'). The direction of kicks (for exam- 
ple, with respect to the spin axis of the neutron star) is 
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more difficult to infer from observations. Recently, several 
observational clues appeared indicating possible NS spin- 
kick alignment. A noticeable spin-kick alignment has been 
inferred from polarization measurements of radio emission of 
pulsa rs (|johnston et al.ll2005l : iRankirJboOTl : Ijohnston et al.l 
i2007i ), as well as from X-ray observations of pulsar wind 
nebu l ae around young pulsars fHclfand . Gotthelf fc HalperrJ 
I2OOII : iKargaltsev. Pavlov fc Garmire, 20061 ). Implications of 
these findi ngs to the formation of d ouble pulsars were dis- 
cussed by IWang. Lai fc Hani l|2006l ). The possibility and 
conditions for such an alignment in the model of the 
kick origin by mul tiple random kicks dur ing NS forma- 
tio n (proposed bv (ISpruit fc PhinnevI Il998l )) were studied 
by IWang. Lai fc Hani (|2007l ). The imphcation of NS kick- 
spin correlat ion to th e plausible birth-kick scenarios was also 
discussed bv iNg fc R omani (2007). 

Here we explore the effect of NS spin - kick correlation 
on the formation and galactic coalescence rate of double 
neutron stars (DNS) which are primary targets for mod- 
ern gravitational wave detectors. We show that the tighter 
alignment, the smaller is the DNS merging rate with re- 
spect to models with random kick orientation. The effect 
is especially important for large kick amplitudes ('^ 400 
km/s). We calculate the spin-orbit misalignment of the com- 
ponents of DNS which can be important for GW data anal- 
ysis. We also considered a scenario in which no (or insignif- 
icant) kick accompanies the formation of a neutron star 
in binary systems from the main-sequence progenitors in 
a restricted mass range (8-11 Mq or so) due to electron- 
capture collapse of O-Ne-Ng degen erate stellar core pro- 
posed bv lPodsiadlowski et al.l l|2004l ) and further elaborated 
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by Ivan den Heuvell (|2004l . l200if ). This hypothesis is phe- 
nomenologically based on the existe nce of long-period Be X- 
ray binaries with low eccentricities ijPfahl et al.ll20o3 ). It is 
consistent with the evolutionary anal ysis of double neutron 
star formation llvan den Heuvei[2007l ) and has been used in 
some population synthe sis studies of DNS, see for example 
iDewi et all (|2005l , (20061 ) . 



2 EFFECT ON THE BINARY NEUTRON 
STAR COALESCENCE RATES 



The effect of NS kick velocity on merging rates of compact 
binaries by means of po pulation synthesis simulations 



was studied earlier (e.g. Il 


Dunov. Postnov & Prokhorov 


llMlt); 


Portegies Zwart & Yungelson 


11998): Belczvnski & Kaloecra I200lf): 


Belczvnski. Kaloffera & BulikI 


[20020") and the inclu- 



sion of the kick into modeling of binary star evolution 
is a prerequisite in all populat i on synthesis simulat ions 
(see IPostnov fc YungelsonI l2006l : iKaloeera et"al l2007l for 
discussion and further references). Clearly, the tight NS 
spin - kick alignment may have important implications to 
the formation and evol ution of b i nary compact stars, as 
was shown previously bv lKalogeral (|2000l ). 

Consider the standard evolutionary scenario leading to 
the f ormation of a binary NS from a mass ive binary sys- 
tem tehattacharva fc van den Heuveil Il99lh. which is also 
discussed in the review iPostnov fc YungelsonI f2006 l. focus- 
ing on the effect of the NS kick velocitjQ. We shall assume 
that the kick velocity vector is confined within a cone which 
is coaxial with the progenitor's rotation axis and charac- 
terized by angle 6 < 7r/2. We shall consider only central 
kicks thus ignoring theoretically fe asible off-cen ter kicks si- 
multaneously affecting the NS spin ([Spruit fc Ph inncv 1998; 
IPostnov fc Prokhorovf [l998l : IWang, Lai fc Hani 1200711 . The 
value of the kick velocity is assumed to obey the Maxel- 
lian distribution f(v) ~ t;^ exp( — (t;/ t)o)^), as suggested by 
pulsar proper motion measurements (iHobbs et al.ll2005h . In 
our analysis we varied the velocity vo from to 400 km/s. 

The rotational axes of both components are assumed 
to be aligned with the orbital angular momentum before 
the primary collapses to form the first NS. The SN explo- 
sion is treated in a standard way as instantaneous loss of 
mass of the exploding star. The effect of the kick on the 
post-explosion binary orbital parameters is treated using the 
energy-momentum conservation in the two point-mass body 
prob l em (see the descriptio n in e.g. iHillj 1 19831 : iKalogeral 
I2OOOI : ICrishchuk et al.ll200j ). The first SN explosion most 
likely occurs when the binary orbit is circular (unless the 
initial binary is very wide so that tidal circularization is in- 
effective), while the second explosion can happen before the 
orbit has been tidally circularized. Possible mass transfer 
phases before the second collapse (such as the common en- 
velope stage and stable mass transfer onto neutron star) are 
assumed to effectively circularize the orbit. In the absence 

^ Here we do not consi der othe r poss ible evolutionary scenario 
for DNS formation, e.g. iBrownl l ll995l) . which allows for hypo- 
thetical hypercritical accretion onto NS in the common envelope. 
These scenarios are studied by means of the popula tion synthesis 
modeling by other authors llDewi et al.ll2005f . I2OO6I) 



of mass transfer the tidal evolution of the orbit eccentric- 
ity is treated according to IZahnI 1119771 ). In our modeling, 
by the time of the second collapse the fraction of eccentric 
binaries which later form DNS attains ~ 10% depending on 
the kick velocity value and direction, as illustrated in Fig. [1] 
It is higher for isotropic kicks and increases with their ab- 
solute values. To treat the explosion in an eccentric binary 
we choose the position of the star in the orbit randomly 
distributed according to Kepler's 2d law. 

We use the population synthesis method 

to calculate the expected coalescence rate of 

DNS (see iIlIpuiiov. Postnov fc Prokhorovl Il997l : 
IPostnov fc YungelsonI l2006f ) and references therein). 
The standard assumptions about binary evolution have 
been made: Salpeter's mass function for the primary's 
mass, dN/dMi ~ M^^'^^, a flat initial mass ratio 
(g = M2/M1 < 1) distribution dN/dq = const, the 
initial semi-major axes distribution in Oepik's form 
dN/d\oga = const. The common envelope phase is 
treated in the standard way based on the energy con- 
servation (|Postnov fc YungelsonI I2OO6I ) with the efficiency 
acE = O.fl The calculations were normalized to the 
galactic star formation rate 3Mq per year, with a binary 
fraction of 50%. The maximum mass of a main sequence 
star which forms a NS in the collapse is set to 30 Mq, 
the maximum mass of a NS is assumed to be 2 Mq. No 
hypercritical accretion onto NS, as assumed to be possible 
in the scenario by Brown (1995), is allowed. We also 
have carefully taken into account rotational evo lution of 
magne ti zed compact stars, as described i n detai l in lLipunovl 
(| 19921 ): iLipunov. Postnov fc Prokhorovl (|l996D . assummg 
no neutron star magnetic field decay. 

The galactic DNS merging rate is shown in Fig. [2] as 
a function of the kick parameter vo and assuming random 
central kicks. The calculations were performed for two as- 
sumptions about kicks - (a) when the formation of NS is 
always accompanied by a kick (we refer to this scenario as 
kick type A) or (b) when the kick is non-zero during the for- 
mation of NS only in binaries starting out from IIA^q , while 
there is no kick velocity at all when a NS is formed in a bi- 
nary system due to the electron-capture core collapse of the 
main-sequence progenitor with mass in the range 8-11Mq 
(kick type B). In the case (a) an almost exponential de- 
crease in the DNS rate coalsecence rate with vo for vo > 100 
km/s is seen. In the case (b) the decrease with vo is less pro- 
nounced, mainly because kickless collapses in binaries are 
more abundant by the assumed Salpeter mass distribution. 

Fig. 13] shows the relative change in the DNS merging 
rate with allowance for the NS spin-kick alignment with dif- 
ferent values of the kick confinement angle 8 for two kick 
models. It is seen that tight alignment (small 6) generally re- 
duces the DNS merging rate, with the effect being especially 
strong for large kick velocity amplitudes. Such a decrease rel- 
ative to calculations with random kicks is clear because the 
NS spin - kick correlation excludes kicks in the binary or- 



^ This important parameter of the evolution of close bina- 
ries is loosely constr a ined, see e.g. the detailed discussion in 
IPostnov fc YungelsonI ||2006| ): however, varying it from 0.1 to 1 
does not change qualitatively the shape of distributions studied 
in the present paper. 
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Figure 1. The distribution of orbital eccentricities before the second collapse in binaries producing DNS for two kick models. Kick 
type A: all NS in binaries receive a kick; kick type B: the NS kick is zero in those binaries where NS is produced from main-sequence 
progenitors with masses 8 — IIMq. 
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Figure 2. Galactic coalescence rate of DNS vs. the kick parameter tio assuming random central kicks. Almost an exponential decay with 
VQ is seen for vq > 100 km/s for kick type A (a), while the decrease in the rate is smaller for kick type B 



bital plane which, if directed opposite to the orbital velocity, 
can additionally bind the post-explosion binary system. 



3 NEUTRON STAR SPIN - ORBIT 
MISALIGNMENT 

There is another observational consequence of the kick in 
DNS systems: the NS spin - orbit misalignment, which can 
be tested by geodetic precession measurements in binary 
pulsars (|Bailes.. .1988). Such a mi salignment is potentiall y 
very interesting for GW studies (|Apostolatos et ai1ll994l ). 
After the first supernova explosion in a binary system (SNl), 
the additional kick imparted to the newborn neutron star 
(NSl) results, in general, in a misalignment between the 



new orbital angular momentum and the NSl (as well as the 
secondary component's) spin vector characterized by some 
angle. For the instantaneous explosion of one of the com- 
ponents treated as point-like masses on a Keplerian orbit, 
this angl e can be calculated analytically, see for example 
iKalogeral l|2000l) . After the second supernova explosion in 
the system (SN2), there are several possibilities for spin- 
orbit misalignment of the compact components. 

1) In close binaries, tidal interactions tend to rapidly 
align the angular momentum vector of the normal star with 
the orbital angular momentum. To spin-up the NS rotation 
to observed ms periods (in binary ms pulsars), a modest 
amount of matter (~ O.IMq) should be accreted by NS. 
This amount is sufficient to align the NS rotation with the 
orbital angular momentum. So if NSl accreted matter before 
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Figure 3. Relative change of DNS merging rate for tlie NS spin-kick correlation as a function of the kick restriction angle d for two kick 
models. 



the second SN explosion, both NSl and the secondary com- 
ponent's spins should be most likely align ed with the orbital 
angula r momentum (see the discussion in I Wang. Lai fc Hani 
lj2006t )l. Note that the NSl spin tends to align with the or- 
bital angular momentum even if NSl does not accrete mat- 
ter but spins-down by the propeller mechanism before the 
second SN explosion, since in that case very strong currents 
must flow through its polar cap and the alignment torque 
can be as strong as during accretion. So in close binaries the 
NSl remains orbit-misaligned prior to the second SN explo- 
sion only in rare occasions where the secondary collapses 
shortly after the first SN in the binary. If both NSl and 
the secondary component were aligned with orbital angular 
momentum prior to the second SN explosion, both neutron 
stars NSl and NS2 will be equally misaligned with the or- 
bital angular momentum (^'i) after SN2 with a kick. 

2) In sufficiently wide binaries, when tidal interactions 
between the components are inefficient, the orientation of 
the NSl spin and the secondary's spin vector may remain 
unchanged until SN2 explosion, after which the orbital an- 
gular momentum vector changes again due to the NS2 kick. 
So in this case we would expect two coaxial NS with spins 
misaligned by angle ^'2 with orbital angular momentum. 
However, such binaries, unless highly eccentric, may be too 
wide to coalesce over the Hubble time. 

We conclude from these considerations that the com- 
ponents of a DNS can have coaxial spins misaligned with 
the orbit by angles ^\ or ^'2 depending on the strength 
of tidal interaction (weak or strong, respectively) acting be- 
tween two SN explosions in the binary system. It is of course 
possible that spins of the components remain misaligned by 
some angle depending on the degree of the spin-orbit inter- 
action of the secondary prior to the collapse. For example, 
NSl spin may conserve its original direction in space, while 
the secondary component before the collapse may have be- 
come aligned with the orbital angular momentum, so in the 
resulting DNS the spin-orbit misalignment angle of the older 
NS will be *2 while that of the younger NS will be *i. In 
this sense, angles ^1 and ^2 should be considered as limiting 
cases. 



Table 1. Mean NS spin-orbit misalignment (in units 7r/2) for 
kick type A. 



i;o. Kick confinement angle 9 (in units 7r/2) 

km/s 0.01 0.1 0.25 0.5 1.0 



50 




0.061 


0.059 


0.057 


0.071 


0.208 






0.307 


0.312 


0.321 


0.305 


0.335 


100 


*i 


0.109 


0.108 


0.113 


0.187 


0.438 




^2 


0.378 


0.380 


0.394 


0.425 


0.629 


200 


*i 


0.190 


0.192 


0.221 


0.311 


0.503 




*2 


0.420 


0.421 


0.456 


0.568 


0.816 


400 


*1 


0.251 


0.257 


0.316 


0.399 


0.547 




*2 


0.454 


0.462 


0.521 


0.671 


0.959 



In our population synthesis simulations we take into ac- 
count the discussed spin alignment effects. In Fig.|4]and[5]we 
show calculated distributions between spins of the compo- 
nents and the orbital angular momentum in coalescing DNS 
systems assuming tidal spin-orbit alignment before the sec- 
ond collapse (angle and the conservation of the original 
components spin direction before SN2 (angle *I'2). 

It is seen that the misalignment angles can be very dif- 
ferent (and even with negative cosines) for random or loosely 
constrained (6 ~ 7r/2) kicks, while a tight spin-kick align- 
men t (6 <g 7r/2) resu lts in much narrow distributions (see 
also iKalogeral l|2000t )'). The mean spm-orbit misalignment 
angles for different values of the kick velocity parameter 
vo and kick models A and B are presented in Table [1] and 
[21 respectively. Fig. H] [S] and Tables [T] [5] also show that the 
difference between vfi and ^2 (and, hence, between spin- 
orbit misalignment of NSl and NS2 in the final pre-merging 
DNS) can be appreciable for small kicks, but tends to be 
less significant for high kick velocities. 

The NS spin-orbit misalignment can be probed by 
studying radio pulsars in binary systems. An extensive anal- 
ysis o f observational data was done by IWang. Lai fc Hani 
(|2006l ). However, high uncertainties in the inferred spin-orbit 
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Figure 4. NS spin-orbit misalignment (cos ^) in coalescing DNS for kick type A with vq = 100 km/s (upper row), vq = 200 km/s (middle 
row) and vq = 400 km/s (bottom row) and different NS spin-kick alignment angles 6. Left panels: the NSl and secondary component's 
spins aligned with the orbital angular momentum prior to the SN2 explosion (the case of close binaries). Right panels: the NSl and 
secondary component's spins aligned with the original binary's orbital angular momentum prior to the SN2 explosion. 
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Figure 5. The same as in Fig. |4] for kick type B. 
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Table 2. Mean NS spin-orbit misalignment 'I' (in units ■7t/2) for 
kick type B. 



^0, 



Kick confinement angle 6 (in units n/2) 



km/s 




0.01 


0.1 


0.25 


0.5 


1.0 


50 




0.039 


0.039 


0.039 


0.054 


0.170 






0.252 


0.254 


0.260 


0.251 


0.270 


100 




0.088 


0.090 


0.092 


0.139 


0.284 




^2 


0.246 


0.250 


0.254 


0.265 


0.391 


200 


*1 


0.167 


0.177 


0.216 


0.282 


0.296 




^2 


0.214 


0.222 


0.260 


0.329 


0.368 


400 


*1 


0.162 


0.191 


0.296 


0.446 


0.510 




*2 


0.176 


0.202 


0.308 


0.458 


0.530 



misalignment angles have not allow firm conclusions to be 
made as yet. 



4 CONCLUSIONS 

We have shown that the spin-velocity correlation observed 
in radio pulsars, suggesting the NS spin-kick velocity align- 
ment, may have important implications to GW studies. 
First, the tight alignment reduces the galactic rate of double 
neutron star coalescences (especially for large kicks 300-400 
km/s) relative to models with random kicks. Second, the 
spin-kick correlation results in a specific distribution of NS 
spin - orbit misalignments. In turn, analysis of the NS spin- 
orbit misalignments inferred from GW signals during DNS 
mergings can be potentially used to put independent bounds 
on the still elusive nature of NS kicks. 
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